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Abstract 

An improved density-matrix expansion is used to calculate the nuclear energy den- 
sity functional from chiral two- and three-nucleon interactions. The two-body interaction 
comprises long-range one- and two-pion exchange contributions and a set of contact terms 
contributing up to fourth power in momenta. In addition we employ the leading order 
chiral three-nucleon interaction with its parameters ce,co and 01^3^4 fixed in calculations 
of nuclear few-body systems. With this input the nuclear energy density functional is 
f— ^ ' derived to first order in the two- and three-nucleon interaction. We find that the strength 

CN . functions Fy{p) and Fso{p) of the surface and spin-orbit terms compare in the relevant 

density range reasonably with results of phenomenological Skyrme forces. However, an 
improved description requires (at least) the treatment of the two-body interaction to sec- 
Q>^ . ond order. This observation is in line with the deficiencies in the nuclear matter equation 

CN i of state E{p) that remain in the Hartree-Fock approximation with low- momentum two- 

and three-nucleon interactions. 

7j : PACS: 12.38.Bx, 21.30.Fe, 21.60.-n, 31.15.Ew 
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\o 1 Introduction 

i^ . The nuclear energy density functional approach is the many-body method of choice in order 

[^ to calculate the properties of medium-mass and heavy nuclei in a systematic manner [H |2]. 

O , Parameterized non-relativistic Skyrme functionals [3l H] as well as relativistic mean-field models 

■ [51 E] have been widely and successfully used for such nuclear structure calculations. A comple- 

r| ■ mentary approach [Zl El El [ID] focuses less on the fitting of experimental data, but attempts to 

constrain the analytical form of the functional and the values of its couplings from many-body 

perturbation theory and the underlying two- and three-nucleon interaction. Switching from 

cd \ conventional hard-core NN-potentials to low-momentum interactions [HI |12] is essential in this 

respect, because the nuclear many-body problem formulated in terms of the latter becomes 

significantly more perturbative. Indeed, second-order perturbative calculations including also 

three-body forces give already a good account of the bulk correlations in infinite nuclear matter 

[T3l [H] and in doubly-magic nuclei p^ . 

In many-body perturbation theory the contributions to the energy are written in terms of 
density-matrices convoluted with the finite-range interaction kernels, and are therefore highly 
non-local in both space and time. In order to make such functionals numerically tractable in 
heavy open-shell nuclei it is desirable to develop simplified approximations for these function- 
als in terms of local densities and currents only. In such a construction the density-matrix 
expansion comes prominently into play as it removes the non-local character of the exchange 
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(Fock) contribution to the energy by mapping it onto a generalized Skyrme functional with 
density-dependent couplings. For some time the prototype for that has been the density-matrix 
expansion of Negele and Vautherin [TB], but recently Gebremariam, Duguet and Bogner [T7] 
have developed an improved version for spin- unsaturated nuclei. They have demonstrated that 
phase-space averaging techniques allow for a consistent expansion of both the spin-independent 
(scalar) part as well as the spin-dependent (vector) part of the density-matrix. The improved 
properties of the new phase-space averaged density-matrix expansion have been extensively 
studied via the Fock energy densities arising from schematic finite-range central, tensor and 
spin-orbit interactions for a large set of semi-magic nuclei (for further details see ref.|17]). 

In order to match with these new developments, the nuclear energy density functional as it 
emerges from chiral pion-nucleon dynamics has been recalculated in ref.[T8]. This calculation 
has treated for isospin-symmetric (i.e. N = Z) nuclear systems the effects from Ivr-exchange, 
iterated Ivr-exchange, and irreducible 27r-exchange with intermediate A-isobar excitations, in- 
cluding Pauli-blocking corrections up to three-loop order. Among other things, it has been found 
that the two- and three-body contributions to the spin-orbit coupling strength Fso{p), as gener- 
ated by 27r-exchange, tend to cancel each other in the relevant density range p ~ 0.08 fni"^, thus 
leaving room for the short-range nuclear spin-orbit interaction. The short-range components of 
the NN-interaction together with the constraints on them provided by the elastic scattering data 
(i.e. NN-phase shifts etc.) have not been considered explicitly in ref.|18j. Furthermore, a simi- 
lar calculation of a microscopically constrained nuclear energy density functional derived from 
the chiral NN-potential at next-to-next-to-leading order (N^LO) has been presented recently by 
Gebremariam, Bogner and Duguet in ref. [I9]. They have proposed that the density-dependent 
couplings associated with the pion-exchange interactions should be added to a standard Skyrme 
functional (with several adjustable parameters). In the sequel it has been demonstrated in 
ref.[2Dj that this new energy density functional gives numerically stable results and that it 
exhibits a small but systematic reduction of the x^-i^ieasure compared to standard Skyrme 
functionals (without any pion-exchange terms). 

The purpose of the present paper is to derive a nuclear energy density functional with 
improved (chiral) two- and three-nucleon interactions. We use for the two-body interaction 
the N^LO chiral NN-potential which reaches at this order the quality of a high-precision NN- 
potential (in reproducing empirical NN-shifts and deuteron properties). The N^LO chiral poten- 
tial consists of long-range one- and two-pion exchange terms and two dozen low-energy constants 
which parameterize the short- distance part of the NN-interaction. The latter contact potential 
written in momentum space gives the most general contribution up to fourth power in momenta. 
In the actual calculation we will use the version N^LOW of the chiral NN-potential developed 
in refs.|2H [22] by lowering the cut-off scale to A = 414 MeV. This value coincides with the 
resolution scale inherent to the universal low-momentum NN-potential Viow-k [HI [I2] to which 
all realistic NN-potentials flow after integrating out effects from momenta above the cut-off 
scale A = 2.1 fm~^. The low-momentum two-body interaction N^LOW is supplemented by the 
leading order chiral three-nucleon interaction with its parameters ce, cd and Ci 3^4 determined in 
calculations of nuclear few-body systems [131 [23]. Our paper is organized as follows. In section 
2 we recall the basic features of the (improved) density-matrix expansion and the nuclear energy 
density functional for isospin-symmetric systems. In section 3 we present the two-body con- 
tributions to the various density-dependent strength functions E{p), Fr{p), Fd{p), -Fso(p) and 
Fj{p), separately for the finite-range pion-exchange and the zero-range contact interactions. In 
section 4, we collect the corresponding analytical expressions for the three-body contributions 
grouped into contact (ce), Ivr-exchange (cd) and 2vr-exchange (01,3,4) terms. Section 5 is devoted 
to the discussion of our numerical results and sections 6 ends with a summary and an outlook. 



2 Density-matrix expansion and energy density functional 



The starting point for the construction of an exphcit nuclear energy density functional is the 
density-matrix as given by a sum over the energy eigenfunctions "^ai"^) representing occupied 
orbitals of the (non-relativistic) many-fermion system. According to Gebremariam, Duguet and 
Bogner [TTj it can be expanded in relative and center-of-mass coordinates, a and r, as follows: 
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with the spherical Bessel function ji{x) = (sinx — a;cosx)/x^. The quantities appearing on the 
right hand side of eq.(l) are: the (local) nucleon density p{f) = 2fc|(r )/37r^ = J2a'^le{^)'^a{r), 

the (local) kinetic energy density T{f) = J2a V\E'a(r ) ■ V\E'Q,(r ) and the (local) spin-orbit density 
J{f) = iJ2a ^l(^)^ ^ V\E'Q,(f^). As shown in section 2 of ref. [12] the Fourier transform of the 
expanded density-matrix eq.(l) with respect to both coordinates a and r defines in momentum 
space a " medium insertion" : 

T{p,q) = l£re-^^^-^[9ikj-\p\) + ^^\kf5'ikf-\p\)-25{kj-\p\) 
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for inhomogeneous many-nucleon systems characterized by the time-reversal-even fields p{r), 
T{f) and J{r). Note that the delta-function S{kf~ \p\) in eq.(2) gives weight to the momentum- 
dependent NN-interactions only in the vicinity of the local Fermi momentum, \p\ = kf{r). 

Up to second order in spatial gradients (i.e. deviations from homogeneity) the energy density 
functional relevant for N = Z even-even nuclei reads: 
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+ Fr{p) 



_2M 4M3 
+ (Vp)2 Fv(p) + Vp ■ JFUp) + P Fj{p) . (3) 

Here, E{p) is the energy per particle of isospin-symmetric nuclear matter evaluated at the local 
nucleon density p{r). The strength function -FV(p) introduces an effective (density- dependent) 
nucleon mass M*{p) and it is related to the single-particle potential U{p, kj) as follows: 
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with p = 2k^r/3iT^. The second equality establishes the equivalent relation to the spin and isospin 

independent p-wave Landau parameter /i(/c/). The strength function -Fv(p) of the (Vp)^ surface 
term has the decomposition [T8] : 



^ , , 1 dFJp) ^ , , 
i^v(p) = 4^ + F.(p) 



(5) 



where -fd(p) comprises all those contributions for which the (Vp)^-factor originated directly 
from the momentum dependence of the interactions in an expansion up to order cp' . Note that 
only the (fixed) nuclear matter piece Qikf — \p\) of the density-matrix expansion goes into the 
derivation of the strength function Fd{p). The second to last term Vp ■ J Fso{p) in eq.(3) is 
responsible for the spin-orbit interaction in nuclei. The associated function Fso{p) measures 
therefore the strength of the nuclear spin-orbit coupling. 
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Figure 1: Finite-range isoscalar and isovector central potentials extracted from N^LOW [22] 



3 Two-body contributions 

In this section we work out the two-body contributions to the various density-dependent strength 
functions which build up the nuclear energy density functional £[p,T,J] written in eq.(3). 
Ideally, one would like to use for this task the universal low- momentum NN-potential Viow-k |12j . 
However, it is generally given in terms of (off-shell) partial wave matrix elements which makes 
its application to the density-matrix expansion rather cumbersome. An explicit representation 
of the momentum space NN-potential in terms of spin- and isospin-operators is much better 
suited for this purpose. For this reason we use (as a substitute for l^iow-k) the chiral NN-potential 
N^LOW developed in refs. [211 122] by lowering the cut-off scale to A = 414 MeV. This value of 
A coincides with the resolution scale inherent to the universal low-momentum NN-potential 
Viow-k- The finite-range part of the N^LOW chiral NNjotential consists of one- and two-pion 
exchange pieces which can be summarized in the form 



■J-pc 

i 



Vr 



M 



NN 



Vciq) + ri-T2 Wciq) + [Vs{q) + n ■ T2 Wsiq)] (Ji ■ <T2 

+ [^r(g) + n ■ ^2 Wriq)] (Ti ■ g <?2 • g 

+ [Vsoiq) + ri-T2 Wsoiq)] ii^i + ^2) -{qxp), 



(6) 



where g denotes the momentum transfer and p the center-of-mass momentum. As usual, ai 2 and 
Ti 2 are the spin- and isospin operators of the two nucleons. A special and simplifying feature 
of V}^j(j- is that all the occurring potentials Vc{q), ■ ■ ■ , Wsoiq) depend only on the momentum 
transfer q and that a quadratic spin-orbit component ~ ai ■ (g x p) c?2 ■ (g x p) is absent. The rel- 
ativistic 1/M^-correction to the 27r-exchange [24J which does (partially) not share this property 
is so small that it can be safely neglected. In order to specify our sign and normalization conven- 
tion, we give also the explicit expression for the Ivr-exchange, Wj^ ^ {q) = ~(5'A/2/7r)^(m^+g^)~"'^, 
with the parameters qa = 1-3, /w = 92.4 MeV and m^ = 138 MeV. 

The solid and dashed lines in Figs. 1,2,3,4 show the finite-range isoscalar and isovector po- 
tentials extracted from the chiral NN-interaction N^LOW [22] in the central, spin-spin, tensor, 
and spin-orbit channel, respectively. In each figure the curves extend up to momentum transfers 
of g = 570 MeV, corresponding to the region g < 2k f within which the interaction gets probed 



■^Note that we associate here the tensor interaction with the operator a\ ■ q S2 ■ q. This operator sphts as 
(g^/3)[S'i2((7) -f o\ ■ 1T2] into the genuine tensor operator S\2{(i) and a spin-spin piece. 
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Figure 2: Finite-range isoscalar and isovector spin-spin potentials extracted from N^LOW [22] 
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Figure 3: Finite-range isoscalar and isovector tensor potentials extracted from N^LOW 
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Figure 4: Finite-range isoscalar and isovector spin-orbit potentials extracted from N^LOW [22] 



for nuclear densities up to p = 0.2 fm^^. One notices in Fig. 3 the large negative values of Wxiq) 
(multiplied with g^/3 in Vj^jy) which result at small momentum transfers from Ivr-exchange. 

In the (first-order) Hartree-Fock approximation the finite-range NN-potential V^vw leads in 
combination with the density matrix-expansion (i.e. by employing the product of two medium 
insertions T{pi,q)T{p2,—q)), to the following two-body contributions to the energy density 
functional £[p,T, J]: 
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setting q = 2xkf. The double-prime in eq.(9) denotes a second derivative and we have given 
the numerical value for F^ip) resulting from the negative curvature of isoscalar central potential 
Vc{q) shown in Fig. 1. One can easily convince oneself that Fj{p) as given in eq.(ll) stays 
finite in the limit fc/ — )■ 0. After expanding the integrand to linear order in g^, the constant 
Vc{0) + 3Wc{0) integrates to zero. 

In addition to the finite-range pieces written in eqs.(7-ll) there are the two-body contri- 
butions from the zero-range contact potential of the chiral NN-interaction N^LOW. The cor- 
responding expression in momentum space includes constant, quadratic, and quartic terms in 
momenta and it can be found in section 2.2 of ref.[25]. The Hartree-Fock contributions from 
the NN-contact potential to the nuclear energy density functional £[p, r, J] read: 
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The 24 low-energy constants Cs,t, Cj and Dj are determined (at the cut-off scale of A 



414 MeV) in fits to empirical NN-phase shifts and deuteron properties 
their values from the pertinent NN-scattering code made available to us 
values of the low-energy constants for the pure contact terms are: Cs '- 
(in units GeV~^), those of the terms quadratic in momenta are: Ci = 
C3 = -29.04, Ci = -524.4, C5 = 717 A, Cq = -42.70, C7 = -1753.6 
those of the terms quartic in momenta are: Di = 0.612, D2 = 25.61, D^ 
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by R. Machleidt. The 
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: 475.8, C2 = 1034.1, 
in units GeV~^), and 
= 19.68, L>4 = -19.49, 
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Figure 5: Three-body diagrams related to the contact (ce) and Ivr-exchange (cd) component 
of the chiral three-nucleon interaction. The short double-hne symbohzes the medium insertion 
r(p, g) for inhomogeneous nuclear matter. 

Ds = 1.287, Dq = 19.02, Dj = 6.565, Dg = -5.429, Dq = 4.226, Diq = -16.02, Du = -1.243, 
Di2 = -0.976, Di3 = -0.998, Du = -7.995, D15 = -0.491 (in units lO^GeV"*^). Let us 
mention that the contributions proportional to Cs,t and Cj in eqs. (12-16) have also been worked 
out in appendix B of ref. [19] and we find agreement with their results. The terms proportional 
to Dj as well as the master formulas eqs. (7-11) for the finite-range contributions are new. Note 
that we do not include an additional regulator function [22] since the NN-interactions are probed 
only at small momenta |pi^2| < ^/ < 285 MeV. 

4 Three-body contributions 

In this section we work out the three-body contributions to the nuclear energy density functional 
S[p,T,J]. We employ the leading order chiral three-nucleon interaction [23] which consists of 
a contact piece (with parameter ce), a Ivr-exchange component (with parameter cd) and a 
27r-exchange component (with parameters Ci, C3 and C4). In order to treat the three-body 
correlations in inhomogeneous nuclear many-body systems we follow ref. [TS] and assume that the 
relevant product of density-matrices can be represented in momentum space in a factorized form 
by r(pi, gi) r(p2, if2) r(p3, — (fi — q2). Such a factorization ansatz respects by construction the 
correct nuclear matter limit, but it involves approximations in comparison to more sophisticated 
treatments outlined in section 4 of ref. [9]. Actually, our approach is similar to the method DME-I 
introduced in ref. [S]. 

4.1 C£;-terni 

We start with the three-body contributions from the contact interaction as represented by the 
left diagram in Fig. 5. With three (inhomogeneous) medium insertions one finds the following 
contribution to the energy per particle: 

which is quadratic in the density p = 2k^/37T'^. Obviously, the contributions to the other strength 
functions F^cLso,j{p) vanish due to the momentum-independence of the contact interaction. For 
the choice of scale A^ = 700 MeV the value ce = —0.625 has been determined in calculations of 
few-nucleon systems [13] (employing in addition Viow-k for the two-body interaction). 






Figure 6: Three-body Hartree and Fock diagrams related to the chiral 27r-exchange three- nucleon 
interaction. 



4.2 CD-term 

Next, we consider the three-body contributions from the Ivr-exchange component of the chiral 
3N-interaction as represented by the right diagram in Fig. 5. Putting in three (inhomogeneous) 
medium insertions one finds the following analytical expressions: 
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with the abbreviation u = kf/rriT^. For the parameter c^ we take the value cn = —2.06 from 
ref.[I3]. Note that there is no contribution to the spin-orbit coupling strength F^oip)- 



4.3 Hartree diagram proportional to ci 3 

We continue with the three-body contributions from the 2'7r-exchange Hartree diagram shown 
in the left part of Fig. 6. Again with three (inhomogeneous) medium insertions one derives the 
following analytical results: 
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which depend only on the two isoscalar couphng constants C\ and C3 with values C\ = —0.76 GeV~^ 
and C3 = — 4.78GeV~^ [23]. Note that the expression for Fso{p) in eq.(25) gives the dominant 
part of the three-body spin-orbit coupling strength suggested originally by Fujita and Miyazawa 
[26] . Their proposed mechanism is based on the excitation of a A(1232)-resonance. In the 
present approach the two-step process vrA^ — ?■ A — ?■ vrA^ is replaced by an equivalent tttiNN 
contact vertex proportional to C3. 



4.4 Fock diagram proportional to ci 3 4 

Finally, there are the three-body contributions from the 27r-exchange Fock diagram shown in 
the right part of Fig. 6. With one single closed nucleon ring this diagram generates (for isospin- 
symmetric nuclear matter) also non-vanishing contributions from the isovector rniNN contact 
vertex proportional to C4. We take consistently the value C4 = 3.96 GeV~^ used in few-body 
calculations by ref. [23]. In the case of the three-body Fock diagram not all of the occurring 
integrals over the three Fermi spheres can be solved analytically. Collecting all the emerging 
pieces, we find the following results for the Fock contributions to the nuclear energy density 
functional £[p,T, J]: 
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Figure 7: Contributions to the energy per particle E{p) of isospin-symmetric nuclear matter. 
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where u = kf/mT^. A good check of all formulas collected in this section is provided by their 
Taylor expansion in kf. Despite the superficial first appearance to the contrary, one can verify 
that the leading term in the /cj-expansion is k'^. In several cases it is even a higher power of kf. 



5 Results and discussion 

In this section we present and discuss our numerical results obtained by summing the series of 
two- and three-body contributions given in sections 3 and 4. The physical input parameters 
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Figure 8: Contributions to the strength function F^-i^p) versus the nuclear density p. 



are: qa = 1-3 (nucleon axial vector coupling constant), /^ = 92.4 MeV (pion decay constant) 
and m^ = 138 MeV (average pion mass). The parameters Cs,t, Cj and Dj corresponding to 
the 24 contact terms in the chiral NN-potential N^LOW [22] have already been listed at the 
end of section 3. For the sake of completeness we quote again the parameters pertinent to the 
chiral three-nucleon interaction: ce = —0.625, cd = —2.06, A^ = 700 MeV, Ci = — 0.76GeV~^, 
C3 = -4.78 GeV-i and d = 3.96 GeV-\ taken from refs. [I3l |23] . 

Fig. 7 shows the contributions to the energy per particle E{p) of isospin-symmetric nuclear 
matter for densities up to p = 0.2 fm~^. The dash-dotted line gives the (attractive) two-body 
contributions and the dashed line the (repulsive) three-body contributions. For comparison we 
have also included the Hartree-Fock contribution to the energy per particle E{p) as obtained 
from the universal low-momentum NN-potential Viow-k [TT] IT2| \TS\ by summing and integrating 
its diagonal (on-shell) partial-wave matrix elements. One observes that our treatment of the NN- 
interaction via the chiral potential N^LOW reproduces these results fairly accurately. The sum 
of the two- and three-body contributions (full line in Fig. 7) shows a first tendency for saturation 
of nuclear matter. However, after inclusion of the kinetic energy E]^i^{p) = 3k'j/10M—3ki/56M^ 
the resulting minimum is still much too shallow. This observation (at the Hartree-Fock level) 
is consistent with refs. fi3\ [H] . An improved description of the nuclear matter equation of state 
E{p) can be achieved when treating the two-body interaction at least to second order. 

Fig. 8 shows the contributions to the strength function Ft(p)- For the two-body part the 
results derived with V^ow-k and the chiral N'^LOW potential he closely together. The three- 
body part (shown by the dashed line in Fig. 8) comes out relatively small. At po = 0.16 fm~^ it 
adds a correction of about 20%. The expression multiplying the kinetic energy density T{f) in 
the nuclear energy density functional eq.(3) has the meaning of a reciprocal density-dependent 
effective nucleon mass: 



M*{p) = M 



3_ 

2M2 



+ 2MF,(p) 



(35) 



It is identical to the so-called "Landau"-mass introduced in Fermi-liquid theory, since it derives 
in the same way from the slope of the single-particle potential t/(p, fc/) at the Fermi surface 



P 



kf. The 



(small) correction term —k'j/2M'^ accounts for the relativistic increase of mass. 



Fig. 9 shows the ratio of effective to free nucleon mass M*[p)/M as a function of the nuclear 
density p. One observes a reduced effective nucleon mass which reaches the value M*(po) ~ 
0.67M at nuclear matter saturation density po = 0.16 fm~^. This is compatible with the range 
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Figure 10: The strength function Fy{p) of the surface term (Vp)^ versus the nuclear density p. 



0.7 < M*{pq)/M < 1 spanned by phenomenological Skyrme forces [3111]. On the other hand 
it has been found recently in ref.[27] that second-order corrections from V|ow-k enhance the 
effective nucleon mass substantially. 

Next, we show in Fig. 10 the strength function Fv(p) of the (Vp)^ surface-term. The pro- 
nounced increase of the two-body contribution (dash-dotted line) at very low densities is caused 
by the Ivr-exchange and has also been observed in other calculations [THl US]- The sizeable 
three-body contribution (dashed line) adds negatively to these initial values such that the total 
result for Fy{p) (shown by the full line in Fig. 10) decreases with increasing density p. For 
comparison we have also included the band (of constant Fv(p)-values) spanned by phenomeno- 
logical Skyrme forces [31 H]. Taking this as a benchmark one sees that our Hartree-Fock result 
is somewhat too small at densities around po/2 = 0.08 fm~^, where the surface energy in nuclei 
gains most of its weight. The calculation of the iterated Ivr-exchange in ref. p,8j suggests that a 
treatment of the low-momentum two-body interaction to second order will further increase the 
values of Fy (p) . 

Of particular interest is the strength function F^oip) of the spin-orbit coupling term Vp • J. 
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The two- and three-body contributions together with their total sum are shown in Fig. 11. The 
two-body part is dominated by the low-energy constant 3C5/8 as indicated also by the weak 
variation of the dash-dotted line with density p. The induced (density-dependent) three-body 
spin-orbit forces add sizeably to this initial value. The major contribution is provided by the 
Hartree term in eq.(23) proportional to C3 = — 4.78GeV~^. With this given value of C3 it is 
considerably larger than the A(1232)-excitation mechanism proposed by Fujita and Miyazawa 
PB] which corresponds to c\t^ = -g\/2/\ ~ -2.9 GeV"^ (with A = 293 MeV the delta-nucleon 
mass splitting). At densities around po/2 = 0.08 fm~^ where the spin-orbit interaction in nuclei 
receives most of its weight, our total Hartree- Fock result overshoots the empirical spin-orbit cou- 
pling strength F^l™^\p) ~ 90 MeVfm^ OH] by about 50%. This requires a compensating effect, 
and indeed it has been found in ref. [18] that the second-order Ivr-exchange tensor force gener- 
ates a spin-orbit coupling of the "wrong-sign". Taking the value FJ^'^^'*''(po/2) — — 35MeVfm^ 
(see Fig. 5 in ref. [IB]) as indicative one can expect that the second-order effects from the low- 
momentum two-nucleon tensor potential will reduce the strength of the spin-orbit coupling 
Fso{p) to the correct amount. 

Finally, we show in Fig. 12 the strength function Fj{p) of the squared spin-orbit density J^ 
in the nuclear energy density functional S[p,T, J]. In contrast to all previous quantities, Fj{p) 
receives only a very small three-body contribution. On the other hand the two-body contribution 
is strongly density- dependent and it reaches quite large values at low densities. This prominent 
feature of Fj{p) has also been observed in previous calculations [HI [19]. Actually, the strong 
density dependence of Fj{p) originates from the dominant Ivr-exchange contribution which 
we reproduce separately in Fig. 12 by the upper dash-dotted line (for an explicit formula for 
Fj ^ (p), see eq.(ll) in ref.[TB]). At this point it should be kept in mind that the J^-term in 
the nuclear energy density functional represents the non-local Fock contributions from tensor 
forces etc. An outstanding Ivr-exchange contribution to the strength function Fj{p) is therefore 
not surprising. 



6 Summary and outlook 

In this work we have used the (improved) density-matrix expansion of ref. [17] to calculate 
the nuclear energy density functional from chiral two- and three-nucleon interactions. We have 
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employed the low-momentum NN-potential N^LOW [22] which is composed of long-range multi- 
pion exchanges and a set of short- distance contact terms. The coefficients of the latter have 
been determined in fits to empirical NN-phase shifts. The leading order chiral three-nucleon 
interaction has been taken with its parameters ce, cd and 01,3^4 fixed in calculations of nuclear 
few-body systems [131 [23] • With this input the nuclear energy density functional £[p, r, J] has 
been derived to first order in many-body perturbation theory, i.e. in the Hartree-Fock approx- 
imation. For the effective nucleon mass M*{p) and the strength functions -Fv(p) and Fso{p) 
of the surface and spin-orbit terms we have found (in the relevant density range) reasonable 
agreement with results of phenomenological Skyrme forces. However, as indicated in particular 
by the nuclear matter equation of state E{p), an improved description of the energy density 
functional requires at least the treatment of the two-nucleon interaction to second order in 
many-body perturbation theory. It is furthermore expected that tensor forces at second-order 
generate an additional "wrong-sign" spin-orbit coupling [ISJ which compensates part of the 
strong three-body contribution to Fso{p)- 

Such a consistent second-order calculation of the nuclear energy density functional repre- 
sents a challenge, since a satisfactory generalization of the density-matrix expansion [K] to 
that situation has not yet been formulated. At second order one must properly account for 
the presence of energy denominators which induce further spatial and temporal non-localities 
and possibly even an orbital dependence of the resulting energy density functional. In a first 
simplified approach one could follow ref.[18j and approximate the energy denominators by the 
kinetic energies of free (on-shell) nucleons, eventually improved by the inclusion of an effective 
nucleon mass M* (p) . In order to keep the second-order calculation manageable and in order to 
see leading effects, one would also restrict the low-momentum two-nucleon interaction to some 
dominant components, such as one-pion exchange with a suitable regularization of its strong 
tensor force at short distances. Studies along this line are underway. 
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